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spacing ~ 0.12 fm, at six sea quark masses m q ranging from m s /6 to m s (where m s is the physical 
strange quark mass). The % t is extracted from the plateau (at large time separation) of the time- 
correlation function of the flavor-singlet pseudoscalar meson (77'X which arises from the finite 
size effect due to fixed topology. In the small m q regime, our result of % t is proportional to m q 
as expected from chiral effective theory. Using the formula Xt = m q^/Nf by Leutwyler-Smilga, 
we obtain the chiral condensate in N f = 2 QCD as Z m (2 GeV) = [252(5) (10)MeV] 3 , in good 
agreement with our previous result obtained in the e -regime. 



The XXV International Symposium on Lattice Field Theory 
July 30 - August 4 2007 
Regensburg, Germany 



* Speaker. 

'Email: twchiu@phys.ntu.edu. tw 



" Physics Department, Center for Theoretical Sciences, and National Center for Theoretical 

Sciences, National Taiwan University, Taipei 10617, Taiwan 
b Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba 305-8571, 

Japan 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribntion-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it/ 



Topological susceptibility in 2-flavor lattice QCD with fixed topology 



T.W. Chiu 



1. Introduction 

In Quantum Chromodynamics (QCD), the topological susceptibility (% t ) is the most crucial 
quantity to measure the topological charge fluctuations of the QCD vacuum, which plays an im- 
portant role in breaking the Ua(1) symmetry. Theoretically, % t is defined as 

Xt = Jd 4 x(p(x)p(0)) (1.1) 

where 

P(x) = ^^£^vXa^[ F tiv{x)F Xa (x)} 

is the topological charge density (oc axial anomaly) expressed in term of the matrix-valued field 
tensor F^ v . 

With mild assumptions, Witten [[I]] and Veneziano [Q] obtained a relationship between the 
topological susceptibility in the quenched approximation and the mass of v\' meson (flavor singlet) 
in full QCD with Nf degenerate flavors, namely, ^(quenched) = f%ml,/(4Nf) where f n = 131 
MeV, the decay constant of pion. This implies that the mass of r]' is essentially due to the axial 
anomaly relating to non-trivial topological charge fluctuations, which can turn out to be nonzero 
even in the chiral limit, unlike those of the (non-singlet) approximate Goldstone bosons. 

Using the Chiral Perturbation Theory (ChPT), Leutwyler and Smilga [^] obtained the follow- 
ing relation in the chiral limit, 

Xt = -^ + #H) d-2) 

where m q is the quark mass, and £ is the chiral condensate. In other words, as m q — » 0, the topo- 
logical susceptibility is suppressed due to internal quark loops. Most importantly, ( |Q| ) provides a 
viable way to extract Z from % t in the chiral limit. 



From (1.1), one obtains 

(Qr) 



* = ^p Qr^Jd 4 x P (x) 



where £2 is the volume of the system, and Q, is the topological charge (which is an integer for QCD). 
Thus, one can obtain % t by counting the number of gauge configurations for each topological sector. 
Obviously, for a set of gauge configurations in the topologically-trivial sector with Q t = 0, it gives 
Xt = 0. However, even for a topologically-trivial gauge configuration, it may possess non-trivial 
topological excitations in sub-volumes. Thus, one can investigate whether there are topological 
excitations within any sub-volumes, and to measure % t using the correlation of the topological 
charges of two sub- volumes [Q]. 

In general, for any topological sector with Q t , using saddle point expansion on the QCD par- 
tition function in a finite volume, it can be shown that ^ (see also |§]) 

lta(pWp(0)>4(|-Z,-^) + ^) (.3) 
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where 

c 4 = -^[(Gf) e =o-3(G?)Lo] 

Thus, one can consider two spatial sub-volumes at time slices t\ and ti, and to measure their 
time-correlation function 



c(h -t 2 ) = (Q(hMt 2 )) = £ (p(*i)p(* 2 )> 



X\,X2 



where the summations run over the spatial volumes at t\ and t2 respectively. Then its plateau at 
large \t\ — tz\ can be used to extract % t , provided that | C4 1 <C 2%}Q.. 

However, for lattice QCD, it is difficult to extract p(x) and Q t unambiguously from the gauge 
link variables, due to their rather strong fluctuations. 

To circumvent this difficulty, one may consider the Atiyah-Singer index theorem [0] 

Q t = n+ = index(^) (1.4) 

where n± is the number of zero modes of the massless Dirac operator & = 7^ (d^ + igA^ ) with ± 
chirality. Since Q is anti-Hermitian and chirally symmetric, its nonzero eigenmodes must come 
in complex conjugate pairs (i.e., Sttp = iX(j) implies £^/50 = —iXj5<p, for X = X* 7^ 0) with zero 
chirality (/ d 4 x^y^ = 0). Thus one can obtain the identity 

n + -n- = Jd\rntr{r 5 ^ +m )-\ X , X )] (1.5) 

by spectral decomposition, where the nonzero modes drop out due to zero chirality. In view of ( |L4] ) 
and (L5), one can regard m tr[y^{& + ni)~ x (x,x)] as topological charge density, to replace p(x) in 
the measurement of % t . 

For lattice QCD, it is well-known that the overlap Dirac operator [j|] in a topological non-trivial 
gauge background possesses exact zero modes (with definite chirality) satisfying the Atiyah-Singer 
index theorem. Writing the massive overlap Dirac operator as 

m\ ( m\ H w (—mo) 



D(m) = U + - + Uo-- y 5 \ 



where H w (—mo) is the standard Hermitian Wilson operator with negative mass —mo (0 < mo < 2), 
then the topological charge density can be defined as 

p l (x)=mti[Y 5 (D c + m)-j c ] 

where (D c +m) _1 is the valence quark propagator with quark mass m, and D c is a chirally sym- 
metric and non-local operator, relating to D(0) by D c = D(0)[1 - D(0)/(2m )] _1 [§]. Note that 
[D c + m) _1 is exponentially-local for sufficiently smooth gauge background and nonzero m. Here 
Pi(x) is justified to be topological charge density, since it can be shown that (see e.g., [|l0|]) 

£pi(x) =mTr[Y 5 (D c + m)- l x ] =n + -n_ (1.6) 
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which is similar to its counterpart in continuum (|Q|), where Tr denotes trace over Dirac, color and 
lattice spaces. 



Now we can replace p(x) with P\(x), and use (1.3) to extract % t for any topological sector. 
However, on a finite lattice, it is contaminated by m n , and any states which can couple to 
(pi(x)pi(O)). A better alternative is to compute the correlator of the flavor-singlet pseudoscalar 
meson 77 ', which behaves as [£|[5|] 



(1.7) 



Then the time-correlation function of r\' (see Fig. [|(a)) is fitted to A +B(e Mt + e M ( T '))toextract 
the constant A = 4^ -% - af 



& - ifn) and provided that |c 4 | < 2x f 2 £l 



c .(t) = 42 





(a) 



o 



disconnected 
connected 

-ti' = disconnect - connect / N f 

A + cosh 
A 




10 15 20 25 

t 



(b) 



Figure 1: (a) A schematic diagram for the time-correlation function of the flavor singlet pseudoscalar me- 
son operator. Each solid line denotes the valence quark propagator (D c + m q )~ l . (b) The time-correlation 
function of the flavor singlet (circles), and its connected (ttiangle down) and disconnected (ttiangle up) 
contributions. Data for m a = 0.025 are shown. 



2. Lattice Setup 

Simulations are carried out for two-flavor (Nf = 2) QCD on a 16 3 x 32 lattice at a lattice spac- 
ing ~ 0.12 fm. For the gluon part, the Iwasaki action is used at j3 = 2.30, together with unphysical 



Wilson fermions and associated twisted-mass ghosts [|11|]. The unphysical degrees of freedom 
generate a factor det[//^(— mo) / (H^—mo) + jj. 2 )] in the partition function (we take mo = 1.6 and 
/I = 0.2) that suppresses the near-zero eigenvalue of H w (—mo) and thus makes the numerical oper- 
ation with the overlap operator substantially faster. Furthermore, since the exact zero eigenvalue is 
forbidden, the global topological change is preserved during the molecular dynamics evolution of 
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the gauge field. Our main runs are performed at Q = 0, while Q = —2 and —4 configurations are 
also generated at one sea quark mass in order to check the consistency as described below. 

For the sea quark mass m q we take six values: 0.015, 0.025, 0.035, 0.050, 0.070, and 0.100 
that cover the mass range m s /6-m s with m s the physical strange quark mass. After discarding 
500 trajectories for thermalization, we accumulate 10,000 trajectories in total for each sea quark 
mass. In the calculation of % t , we take one configuration every 20 trajectories, thus we have 500 
configurations for each m q . For each configuration, 50 pairs of lowest-lying eigenmodes of the 
overlap-Dirac operator D(0) are calculated using the implicitly restarted Lanczos algorithm and 
stored for the later use. 



3. Results 

For the connected diagram (see Fig. 0(a)), the pion correlator is computed using the conju- 
gate gradient algorithm with a low-mode preconditioning. Low-modes are also used for averaging 
over source points [O], which significantly improves the statistical signal. For the disconnected 
diagram, the quark propagator is represented by the eigenmode decomposition and approximated 
by the 50 conjugate pairs of the low-lying eigenmodes. The quark propagator is then obtained for 
any source point without extra computational cost, and the disconnected loops can be calculated 
with an average over the source point. The truncation is motivated by the expectation that the long 
distance correlation is dominated by the low-lying fermion modes; its validity has to be checked 
numerically (see below). 

In Fig. |l|(b), we plot C^{t) together with those of connected and disconnected parts for m q = 
0.025. The curve is a fit to a function A + B(e~ Mt + e~ M ^ T ~'^ ) with data for C n >{t) in the range 
t € [4,28]. The horizontal line is a fitted constant A = 1.70(13) x 10~ 3 , where the error is estimated 
using the jackknife method with bin size of 20 configurations, with which the statistical error 
saturates. Assuming \c$\ <C 2;^ 2 £2, we obtain a 4 Xt = 3.40(27) x 10~ 5 at m q = 0.025. 

Since the disconnected diagram is computed with only 50 pairs of low-lying eigenmodes, we 
have to check whether they suffice to saturate C^{t). For the time range [4,28] used for fitting, as 
the number of eigenmodes is increased from 10 to 30, the change of correlator \8C r] i\/C r] i is ~ 3%, 
while from 30 to 50, it is only ~ 0.3%, which is less than 8% of the statistical error. Thus C^> is 
well saturated with 50 eigenmodes. This also holds for all six sea quark masses. 

In Fig. ||, we plot the topological susceptibility % t a 4 as a function of the sea quark mass 
m q a, The statistical precision is good enough to find a clear dependence on the sea quark mass. 
For the smallest three quark masses, 0.015, 0.025, and 0.035, the data are well fitted by a lin- 
ear function F + Gm with the intercept F = 0.0(1) x 10~ 5 and the slope G = 0.00133(5). Evi- 
dently, the intercept is consistent with zero, in agreement with the #PT expectation ( |1.2[ ). Equat- 
ing the slope to a 3 H/Nf, we obtain a 3 L = 0.0027(1). In order to convert £ to that in the MS 
scheme, we calculate the renormalization factor Z ; ^ s (2GeV) using the non-perturbative renor- 
malization technique through the RI/MOM scheme Q. Our result is Z,^ s (2GeV) = 0.742(12) 
Q. With a- 1 = 1670(20) (20) MeV determined with r = 0.49 fm the value of I is tran- 
scribed to £ MS (2GeV) = [252(5)(10)MeV] 3 , which is in good agreement with our previous result 
[251 (7) (11) MeV] 3 |jT§] obtained in the e-regime from the low-lying eigenvalues. The errors rep- 
resent a combined statistical error (a -1 and Z ; ^ s ) and the systematic error estimated from the higher 
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Figure 2: Topological susceptibility % t a A versus sea quark mass m„a. 



order effects (e.g., C4 term), respectively. Since the calculation is done at a single lattice spacing, the 
discretization error cannot be quantified reliably, but we do not expect much larger error because 
our lattice action is free from 0(a) discretization effects. Our results of topological susceptibility 
are listed in Table |], for six sea quark masses respectively. 



m q r 


0.0616(8) 


0.1016(15) 


0.1412(17) 


0.1982(22) 


0.2742(28) 


0.3852(36) 


xA 


5.59(47) x 10~ 3 


9.26(1.29) x 10~ 3 


0.0123(22) 


0.0187(34) 


0.0263(33) 


0.0199(33) 



Table 1: The values of topological susceptibility % t rQ extracted in the 2 = sector, for six sea quark masses 
respectively. 

In principle, Xt in ( |1.3| ) is universal for any topological sector. We check the universality of 
Xt as follows. At sea quark mass m q = 0.050, we generate 250 configurations with Q = —2 and 
—4 respectively in addition to the main run at Q = 0. Then we extract Xt from the time-correlation 
function of r\', similar to the Q = case. Our results for a'Xt are: {7.4(1. 3), 6.4(2.1), 5.9(1. 8)} x 
10~ 5 for Q = {0, —2, —4} respectively. Evidently, Xt extracted from different topological sectors 
are consistent with each other within the statistical error. 

Finally, we come to the assumption \c$ \ <C 2#, 2 n in extracting Xt vm ( p"^ )- With the formulas 
derived in j5|], we can obtain an estimate of the upper bound of \c$\/ (2x?Q) by measuring the 
two-point correlator (pi(jfi)pi(xi)) and the four-point correlator {pi(x\)p\(xz)p\(x-i)p\(x4)) . Our 
preliminary result is \a{\/ (2x?&) < 0.1, for all six sea quark masses. Details of this calculation 
will be presented elsewhere. We note that the upper bound \cn\/ (2x}&) < 0.1 is also consistent 
with the ratio \ct\j%t ~ 0.3 obtained in the quenched approximation J17Q. 



4. Conclusion 



In this paper, we have determined the topological susceptibility Xt m two-flavor QCD from a 
lattice calculation of two-point correlators at a fixed global topological charge Q t = 0. The expected 
sea quark mass dependence of Xt from #PT is clearly observed with the good statistical precision 
we achieved, in contrast to the previous unquenched lattice calculations. Our result indicates that 
the topologically non-trivial excitations {e.g., instanton and anti-instanton pairs) are in fact locally 
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active in the QCD vacuum, even when the global topological charge is kept fixed. The information 
of these topological excitations is carried by low-lying fermion eigenmodes if the exact chiral 
symmetry is preserved on the lattice. This work demonstrates that Monte Carlo simulation of 
lattice QCD with fixed topology is a viable approach, to be pursued when the topology change 
hardly occurs near the continuum limit even with chirally non-symmetric lattice fermions. The 
artifacts due to the fixed topology in a finite volume can be removed to obtain the physics results 
in the 6 vacuum, provided that % t has been determined in the first place [0,0], as has been done in 
this work. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 07-16), and also in part on NEC SX-8 at YITP (Kyoto U), NEC SX-8 at 
RCNP (Osaka U), and IBM/HP clusters at NCHC and NTU-CC (Taiwan). This work is supported 
in part by the Grant-in-Aid of the Japanese Ministry of Education (No. 13135204, 15540251, 
17740171, 18034011, 18340075, 18740167, 18840045, 19540286, and 19740160) and the Na- 
tional Science Council of Taiwan (No. NSC95-2112-M002-005, NSC95-2112-M001-072). 
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